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Q ■ Abstract 

^ ' The Hubbard model of interacting electrons, like the Ising model of spin- 

lO ■ spin interactions, is the simplest possible model displaying many "real world" 

features, but it is much more difficult to analyze qualitatively than the Ising 
(N . model. After a third of a century of research, we are still not sure about many 

of its basic properties. This mini-review will explore what is known rigorously 



> 



, about the model and it will attempt to describe some open problems that are 

possibly within the range of rigorous mathematical analysis. 



0. Introduction 

The Hubbard model is to the problem of electron correlations as the Ising model is 
I to the problem of spin-spin interactions; it is the simplest possible model displaying many 
Q ' "real world" features. It is, however, much more difficult to analyze qualitatively than the 

J-^ ■ Ising model. After a third of a century of research, we are still not sure about many of its 
basic features. It is believed, for example, to have something to do with high temperature 
^ ■ superconductivity, and it would be nice to validate this. 

■ One of the most intriguing questions concerns magnetism — particularly the magnetic 

properties of the ground state. The familiar models of interacting spins, such as the Ising 
or Heisenberg models, posit an underlying spin Hamiltonian that ultimately comes from 
a model of itinerant electrons. The mystery here is that neither the electronic kinetic 
nor the electronic potential energies favor ferromagnetism, but together they sometimes 
do so. Ferromagnetic spin models, for example, are popular, but the truth is that the 
antiferromagnetic preference of the kinetic energy usually seems to dominate. In fact 
the only known examples of saturated itinerant electron ferromagnetism come from the 
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Hubbard model in a special limit (infinite repulsion and one hole) devised by Nagaoka [NY] 
or from the Hubbard model with finite repulsion but on special lattices devised by Mielke 
and Tasaki [Mil, TH2, MT] for which the kinetic energy spectrum is macroscopically 
degenerate. No doubt, further study of the Hubbard model will eventually lead to a better 
understanding of the way in which the Pauli exclusion principle leads to magnetism. 

These notes will explore what is known rigorously about the model and it will attempt 
to describe some open problems that are possibly within the range of rigorous mathematical 
analysis. The notes are not a complete review in the sense that every contribution is 
covered, but an attempt has been made to mention, at least, most topics that can be 
treated rigorously. I apologize to those authors whose works have been omitted and I hope 
this will be attributed to my ignorance rather than intent. Given more space I would have 
liked to discuss the closely related Falicov-Kimball model, for which many rigorous results 
are known, but which is not a Hubbard model because it does not have SU (2) symmetry. 
Another interesting chapter would be the infinite dimensional Hubbard model — pioneered 
by Metzner and VoUhardt [MV] — and which opens exciting mathematical and physical 
avenues. It has led to a large literature, but much remains to be added in the way of 
mathematical rigor. 

I thank Peter Eckle, Walter Metzner and Hal Tasaki for critically reading this manuscript 
and I thank Balint Toth for considerable help with the bibliography. 

1. Definition of the Model 

Only the original short-range Hubbard model will be considered here. Like the Ising 
model, the Hubbard model is defined on a graph, i.e., a collection of vertices or sites 
(denoted by A and whose number is |A|) and (unordered) edges or bonds connecting 
certain distinct pairs of vertices. The word graph instead of lattice is used to avoid any 
possible implication of translation invariance because most of the results stated here do 
not depend on such invariance. One is given a hopping matrix T, with elements txy, 
with X and y G A, and we assume, as a convention, that txy — ii x,y are not connected 
by an edge. Note that txx = 0. txy might be complex, signaling the presence of a magnetic 
field, the line integral of whose vector potential from x to y (thought of now as points in 
R^) is eLTg{txy). However, T is always self adjoint, txy = t*^, with * denoting complex 
conjugate. 

The bipartite graphs form an important sub-class; here A = AU B, with A and B 
disjoint, and such that there is no edge between x,y if x E A and y G A or if x E B and 
y E B. The square lattice is bipartite, the triangular is not. 

Electrons, i.e., spin ^ fermions, move on A with kinetic energy given in second quan- 
tized form by = K-^ + Ki with 




(1.1) 



x,yeA 
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Here cr = ±1 denote the two spin states t and J, while c^^ is the creation operator for an 
electron at x with spin cr. We have c^o-St ~'~ ^yr^lia ~ ^xy^ar and c^aCyr + CyrCxa = 0. 
The number operator is defined by Uxa = c^ct^xct eigenvalues and 1. The total 

number of each spin species 

xeA 

is a conserved quantity. The total particle number is N = N-^ + which satisfies < 
N < 2|A|. The half-filled band, N = |A| is especially important and especially amenable 
to analysis. 

In the physics literature it is often assumed that t^y = constant = t > on all edges of 
A, in which case T is a discrete version of the Laplacian, but without the diagonal terms. 
This assumption will not generally be made here. 

At each site x there is also given a number 17^ which governs the on-site electron- 
electron interaction at x. We usually assume all > (repulsive case) or all Ux < 
(attractive case). The total potential energy is''' 

W=J2Ux{nx^-k){nxi-^). (1.2) 
xeA 

The total Hubbard Hamiltonian of our system is then 

H = K + W, (1.3) 

and most of the results discussed here are about the ground state of this H. In fact, 
many of the results do not require a point interaction as in (1.2); instead, terms like 
Uxyiux-] — ^)(^yi ~ with the matrix Uxy being positive semidefinite, are also allowed. 

The interaction in (1.2) includes a one-body term — i Ux{nx^ + rixi) plus a trivial 
constant term ^ Ux- If Ux is independent of x, as is normally assumed in almost all 
papers on the subject, then this one-body term is trivially a constant proportional to the 
fixed particle number. It is interesting to consider Ux ^ constant, and the one-body term 
is included in (1.2) in order to be able to exploit hole-particle symmetry: The unitary 
transformation that maps c^^ — > c^^ and cl^ — > c^^ (but c^-f Cxi) is the hole-particle 
transformation on the down-spins and it maps (uxi — ^) into — (n^j, ~ ^)- Thus, this 
transformation maps the repulsive W into the attractive W and vice versa, but it changes 
the down-spin number from to |A| — N^. 

This hole-particle transformation may or may not map Ki into Ki. It does map Ki 
into Ki, defined by txy = ~t*^. 

^ The usual formulation is n-^n^ instead of our (n^ — ^)(nj — |). If Ux varies with x 
the two formulations are obviously inequivalent. The formulation here is not without a 
physical foundation because a neutral atom (with n^^ -|- nj^ = 1) is locally the most stable 
configuration and adding or removing an electron produces a net local charge that raises 
the energy roughly equally. 
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In the special, but important case that T is real and A is bipartite we can make a 
further unitary transformation that will take into and Uxi into Uxi- This unitary 
transformation maps Cxy into {— 1)^0x1 (^-^id, of course, Cx] into c^-f), where (—1)^ denotes 
the function on the vertices of A which is +1 for x E A and —1 for x & B. Effectively, txy 
is mapped into (— 1)^(— l)^^^;^ = —txy for the down-spins. Henceforth, the hole-particle 
transformation on bipartite lattice is always meant to include this additional unitary (— 1)^, 
so that Ki is mapped into . 

Thus, in the special real, bipartite case we can analyze the repulsive case by analyzing 
the attractive case. But we have to remember that ^ |A| — A'"j^ so that the particle 
number for one is related to the magnetization of the other. This is a non-trivial distinction 
and is similar to the relation between the ferromagnetic and antiferromagnetic Ising models 
on a bipartite graph. 

The hole-particle transformation also induces a conservation law in the real, bipartite 
case that has no classical analogue. Like any electron system this model has an SU{2), 
i.e., angular momentum, invariance. The generators are 

j3 = 1 (iV^ - NO, J+^Yl 4^.^ J- = (J^)^- (1-4) 

xeA 

However, the Hamiltonian i7 in (1.3) is the unitary transform (described above) of another 
H' (with Ux replaced by —Ux) which also has an SU{2) symmetry with generators given 
by (1.4) in the transformed basis. By transforming back we obtain a different set of SU (2) 
generators 

J^=l{N^ + N^- |A|), J+ = E(-ir4T4i, J- = (J'-y (1-5) 

that also commute with H and with the J's in (1.4). We call this the pseudospin. This 
special model thus has an SU (2) x SU (2) symmetry group, ^ but it must be emphasized that 
real T is essential (zero magnetic field). These operators change the quantum numbers 
and iVj. The operators change by one unit while change by one unit. The four 
operators together thus permit us to move around inside rectangles in (N-^, Ni) space whose 
four vertices have the form (n,m), (m, n), (|A| — n, |A| — m), (|A| — m, |A| — n). This means 
that to each eigenstate for = n, = m there is a corresponding eigenstate with the 
same energy at each point in the rectangle (there are also additional eigenstates, of course). 
Thus, the usual operators permit us to infer all the states of our Hamiltonian from 
knowledge of the states on the line N^^ — = or the line — = 1. The pseudospin 

§ Actually, the group is 50(4) = SU (2) x SU (2)/Z2, as pointed out by Yang and Zhang 
[YZ]. The reason is that the two operators w = — 1 in each SU{2) corresponds to only 
one operator on our Hilbert space, i.e., w ®\ = \ ®w. This reduction to SO{4) coincides 
with the observation that in every state the spin and pseudospin are either both integral 
or both half-integral. 
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operators J permit us to infer everything from knowledge of the hne + Ni = |A| 
(half-filled band) or the line A?"-!- + A?"! = |A| + 1. But we repeat that this property of 
holds only for the real, bipartite case. 

We note, for future use, that the hole-particle transformation {without (—1)^) applied 
to both spins (namely c^^ ^ cj.^). does not preserve J. Indeed, — > —J^,J^ — > — J^, 

but J2 ^ J2. 

The Hubbard model describes — in the simplest possible fashion — an interacting 
fermion system. It can be viewed this way, as a toy model, or it can be viewed, a bit 
more realistically in the repulsive case, as a serious model of 7r-electrons hopping between 
localized Wannier orbitals in some molecule such as benzene (with |A| =6); the half-filled 
band, N — |A|, is then especially important because it corresponds to neutrality. The 
ultra-short range interaction is supposed to mimic a highly screened Coulomb potential. 
Prom the latter viewpoint it was known first in the chemistry literature as the Pariser- 
Parr [PP]-Pople [PJ] model; molecules having a bipartite structure are called "alternant 
molecules". It was a decade later that Hubbard [HJ], Gutzwiller [GMC] and Kanamori 
[KJ] realized its importance for bulk matter. 

2. One-Dimensional Exactly Solvable Model 

In 1968 the ground state of H was solved for the translation invariant one-dimensional 
ring by Lieb and Wu [LW] using the extension of the "Bethe ansatz" technique [BH, LLl] 
to fermions [MJ, FLl, YC, GM]. Shortly thereafter, Ovchinnikov [OA] used these results to 
calculate the elementary excitation spectrum at half filling. Recently, Essler and Korepin 
[ES] obtained a new and illuminating derivation. Coll [CC] extended Ovchinnikov's results 
to arbitrary filling, as did Woynarovich [WF]. Takahashi [TM] evaluated the magnetism 
curve at half filling and Shiba [SH] evaluated the magnetic susceptibility for all filling. 
Koma [KO] has formulated a sequence of approximations to the positive temperature free 
energy and correlation length (obtained by using the Trotter product formula for the 
partition function) which converge to the exact answer and such that each approximation 
can be calculated using the Bethe ansatz without any assumptions, such as the "string 
hypothesis" . 

These results were obtained in the thermodynamic limit, in which sums could be 
replaced by integrals. Woynarovich and Eckle [WE] evaluated the asymptotics of finite 
size effects on the ground state energy. For small chains the general "Bethe ansatz" 
solution, while correct, is too complicated for numerical evaluation. Heilmann and Lieb 
[HL] undertook to evaluate all the energy levels for all U > for the benzene molecule 
{N = 6, |A| = 6). To our surprise we found many instances of both level crossings and of 
permanent degeneracy — as a function of U — and which were not accounted for on the 
basis of the known invariance groups (spin, pseudospin and symmetries of the hexagon). 
This means that the system has non-abelian symmetry groups and these are dependent 
on U (i.e., the group operations commute with H but not with K and W separately; a 
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well known example of this phenomenon, for the hydrogen atom, is the Runge-Lenz vector 
whose definition depends on the value of the electron's charge). 

What are these invariants? Sixteen years later Shastry [SB] (unaware of [HL]) found 
many invariants and also a two-dimensional classical statistical mechanics vertex model 
whose transfer matrix commutes with our H. In fact he found a whole commuting family 
of such transfer matrices — which means that the one-dimensional model can be called 
"integrable" . It is not clear whether all invariants of H are of Shastry's form. A few 
years later, Grosse [GH] (motivated by [HL]) published another derivation of some of the 
invariants. 

Another interesting question is whether all the eigenstates of H arc of the "Bethe 
ansatz" form. This was answered in the negative by Essler, Korepin and Schoutens [EKS] 
who went on to demonstrate that when the ring is bipartite (i.e., even length) the SU{2) x 
SU{2) generators and </+ save the day. The lowest weight states (the ones that are 
annihilated by J~ and J~) are claimed to be all of the "Bethe ansatz" type. All the 
remaining states are then obtained by application of J"*", and J"*". 

Many of the "Bethe Ansatz" results rely on a "string" hypothesis. Moreover, some 
assertions in [EKS], notably the linear independence of the solutions, rely on an appeal to 
some properties of Shastry's invariants that have not been verified. It is desirable to put 
these matters on a more rigorous basis. 

The literature about this one-dimensional solution is vast and the above remarks do 
not reflect everything that is known about the subject. 

3. Magnetism 

A. One-Dimension. 

It is convenient, now, to take an open chain instead of a ring. Then, by a general 
theorem of Lieb and Mattis [LMl] , whose proof also works for lattice systems with nearest 
neighbors hopping and a completely arbitrary many-body potential, (S), the ground 
state energy of H, as a function of total spin S < and for N particles, satisfies 

E^'ilN) > E^(|Ar - 1) > E'^i^N - 2) > . . . > £;^(0) or £;^(|). (3.1) 

The numbers are now totally arbitrary and not necessarily of one sign. [Note: The 
original theorem [LMl] uses a Perron-Frobenius positivity argument, which a-priori re- 
quires txy > 0. However, any T can be reduced to this form in one-dimension by 
the simple unitary gauge transformation c^a exp[— i6'(x)]ca;o-, which carries t^y into 
txy '•= ^a;j/ exp[z6'(a;) — iO{y)], where 9{x) is an arbitrary real number, selected so that 
txy > for all X and y.] Since we can then take T to be real, our H is then the (down- 
spin) hole-particle transform of some other H, in which is replaced by —Ux- Assuming 
to be even, (3.1) says that the ground state of H has N-^ = = N/2. The transformed 
H system has N-^ = N^, Ni = |A| — N^, whence = ^{N-f — Ni) is the predetermined 
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number ^{N — |A|). The H system also satisfies (3.1) since it is one-dimensional, and the 
spin of the ground state of the H system is therefore 

||iV- |A||. 

This number, S, is the pseudospin of our ground state — a curious result whose physical 
significance is not entirely obvious. (Note the logic here. Both systems, H and H must be 
in their respective ground states, consistent with the given conditions on each; for H it is 
N^+Ni = N and for H it is - = N - |A|.) 

This theorem was extended many years later by Aizenman and Lieb [AL] to posi- 
tive temperatures. The main theorem, applicable to any many-body potential in a one- 
dimensional system, expresses the fact that the free energy is a monotone increasing func- 
tion of the spin. This is done in terms of total spin S or 3-component s, as follows; neither 
inequality implies the other. 

Classify the eigenstates by the 3-component of spin, = s = ^{N-^ — N^), and by 
5", the total spin angular momentum (recall = S{S + 1)). For a given and inverse 
temperature /3, let Zs{s) = Tr^e"^^ be the partition function in which only states of a 
given s value are included. Likewise, let Zj{S) = Tr^e"^^ be the partition function for a 
given S. We have the relation Zz{S) — ^3(6'+ 1) = (25*+ l)^^Zj(S'), which is obvious from 
the theory of angular momentum. Correspondingly, we have the combinatorial quantities 
Yz{s) := (jv^J and yj(,S) := (2,S+l)[y3(;S)-y3(;S+l)], which are essentially the partition 
functions of free particles with iif = 0, and which serve to normalize the Z's. The two 
theorems are then that 

Ysis) YAS) 

are both strictly monotone decreasing functions of their arguments {s or S). A corollary of 
this is that the magnetization is less than its value in the atomic limit (or what I prefer 
to call the pure paramagnetic value). I.e., for all /? and magnetic field h 

M{/3,h)^^^lnZ{h) <Ntaiih{ph). (3.2) 

LJ CI/ t 

N/2 

Here Z{h) = ^ Z3{s) exp[2/3s/i] is the total partition function (recall that the ^r-factor 

s=-N/2 

of an electron is 2). 

It is noteworthy that all this holds for completely arbitrary, Hermitian txy and (real) 
Ux- Indeed, it holds even if we add an arbitrary real one body potential 

V = y^^Vxirixt + na;i). 

An amusing fact about one-dimension concerns the U = +00 case. For an open chain 
the ground state is highly degenerate — indeed, it can have any value of S. For a closed 
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chain the situation is quite different and the S of the ground state depends on the sign of 
the txy's and whether is even or odd; S can be N/2 in some cases [AL]. The situation is 
discussed in detail by Mielke [Mi2] who finds that the "average" 5" is \/iV- Some interesting 
facts about the closed one-dimensional chain with [/ = 0, but with a magnetic field, can 
be found in [LL2]. The U case is discussed by Fujimoto and Kawakami [FKl]. 

B. Half-Filled Band. 

The hole-particle symmetry notwithstanding, the repulsive case (all Ux > 0) and the 
attractive case (all Ux < 0) are quite different, even for a bipartite lattice. The physical 
spin of one is the pseudospin of the other. 

In the limit U — > -|-oo (by which I mean all Ux —>■ +oo), the energy and wave functions 
have nice limits. The electrons become hard-core particles. When N = |A| we just have 
one electron per site and, since motion is impossible, each electron can be, independently, 
spin-up or spin-down. This gives us the atomic limit whose partition function is YJ^ or Yj 
given above. First order perturbation theory in t/U vanishes, but in second order wc have 
to diagonalize our H among all the 2l^l degenerate ground states just described. This 
yields an effective Hamiltonian H' that can be written in terms of the three Pauli spin 
operators S at each site. It turns out [AP] that for any graph 

H' = Jxyi^x ■ — |) (3-3) 

x,yeA 

with Jxy — \txy\'^{Ux^ + Uy^). This is the spin 1/2 antiferromagnetic Heisenberg Hamil- 
tonian and it is known [LEI, LM2] that its ground state has total spin 

,S= - (3.4) 

on a bipartite graph. It is also known [DLS, KLS] that in the translation invariant 
case (i.e., our graph is a D-dimensional hypercube with periodic boundary conditions and 
txy = constant = t and Uxy ~ constant = U) there is long range order when D > 3. 

The obvious question is whether the results just stated (i.e., (3.4) and the long range 
order) hold non-perturbatively in the repulsive case. One would also guess that in the 
attractive case the total spin should be zero in the ground state because when Ux = —oo 
for all X the ground state consists simply of bound pairs of electrons sitting on selected 
sites. These questions about the spin are answered in the following [LEI]. 

Theorem 1: Assume txy is real for all x,y G A. If Ux < for all x, the ground 
state on any connected graph is unique and has spin S' = for any even electron number, 
not just AT = |A|. If Ux > Q for all x, if A is connected and bipartite and if N = |A| 
is even, the ground state is unique (except for the {2S + l)-fold degeneracy) and has spin 
S=l\\A\ - \B\\. 

We can easily have |^| — \B\ of the order of |A| itself. As an example, take a square 
lattice and add a site at the center of each bond of this square lattice. The original sites 
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are then B sites and the new sites are A sites. Then |^| = 2\B\ and the ground state 
has a bulk magnetization per site of 1/3. This is really more like ferrimagnetism than 
ferromagnetism but, in any case, it is one of the few examples known in which the system 
has a bulk magnetization without an external magnetic field. 

There is an interesting corollary of this theorem if A is bipartite and if l^l > \B\. 
Suppose we ask for the absolute minimum energy, without fixing A^. Starting with Ux > 
we find that an optimum N is N = |A| (by the SU{2) x SU{2) symmetry mentioned in 
Sec. 1). But J3 = (Ar^-Ar|)/2 can be anywhere in the interval (|5| - 1^|)/2 to (|^|-|S|)/2 
since S = {\A\ — \B\)/2. Now using hole-particle symmetry to obtain the < model, we 
find that the optimum A^ is any integer satisfying 2\B\ < N < 2\A\. (E.g., starting with 
N]- = \A\,Ni = \B\, the transformed values are A"-|- = A"j = |A| — \B\, which yields 
A^ = 21^41. ) Thus, there can be a large degeneracy in the attractive case! (For U — this 
is easily seen from the remark that T has at least \A\ — \B\ zero eigenvalues.) 

This theorem was extended to positive temperature by Kubo and Kishi [KK] who 
found upper bounds on certain two-point functions. They discuss only the translation 
invariant case on a hypercubic lattice with Ux = U = constant, but their method easily 
extends to the general case. For U < they bound the spin susceptibility at wave vector 
q by 

X, < (3.5) 

for all temperatures and all filling fraction A^/|A| (more precisely, they use the grand 
canonical ensemble and prove (3.5) for all chemical potentials /i). This result precludes 
magnetic long range order. In the repulsive case U > 0, and with /j, adjusted to the 
half- filled band case N = |A|, namely = 0, they bound the charge susceptibility as 

{5nq,5n-g)<{PU)-^ (3.6) 



and the pairing susceptibility as 

{p„p-,) < {pU)-\ (3.7) 

Here {A,B) is the Duhamel two-point function {A, B) = TrfA+e^^^Se^^"*)^^] dt and 
5nq = Uq — (riq) aud Px = cl^Cxi and ^ denotes spatial Fourier transform. Charge long 
range order is precluded by (3.6) while Cooper pairing is precluded by (3.7). 

Although I am restricting this review to the Hubbard model, I cannot resist the temp- 
tation to mention that Theorem 1 has recently been extended [FL2] to another model — the 
Holstein model — in which electrons interact with a quantized phonon field instead of with 
each other. Again, the finite system ground state for any even number of electrons is 
unique and has zero spin. The method of proof of Theorem 1 has also been used by Ueda, 
Tsunetsugu and Sigrist [UTS] to show that the periodic Anderson model at half-filling has 
a singlet ground state. 
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This close connection between the half-filled band, repulsive Hubbard model and the 
antiferromagnetic Heisenberg model points to the first of our spin problems. The antifer- 
romagnetic Heisenberg model on a hypercubic lattice, in the thermodynamic limit, has no 
long range order (LRO) at positive temperature in dimensions D = 1 or 2. This is a con- 
sequence of the Hohenberg-Mermin- Wagner theorem. The same is true for the Hubbard 
model as first shown by Walker and Ruijgrok [WR], then by Ghosh [GD]. Later, Koma 
and Tasaki [KT] proved it by a different method — that of McBryan and Spencer [MS]. As 
far as the ground state is concerned, the Heisenberg model has LRO for D = 2 and spin 1 
or more per site [KLS] (i.e. |Sj;|^ = S{S + 1) and S > 1) and it is believed to have LRO 
also when S = 1/2, which is the case of interest for us. The D = 1 case is believed to 
have no LRO in the ground state. For D > 3 there is LRO in the ground state for all 
> 1/2 [KLS] and at positive temperature [DLS, KLS] for S' > 1 (presumably, also for all 
S > 1/2). The obvious conjecture is the following. 

Problem 1: Prove that there is antiferromagnetic LRO in the half-filled band, re- 
pulsive Hubbard model on the hypercubic lattice (with t^y = 1 for \x — y\ = 1 and Ux = 
positive, finite constant) in the ground state when D = 2 and for positive temperature when 
D >3. For which values of N will the attractive model have LRO? 

The Falicov-Kimball [FK2] model poses an analogous problem that can be solved af- 
firmatively. In this model, one kinetic energy term, Ki, is omitted from the Hamiltonian. 
Although the down spins are not dynamic their locations are left arbitrary. After "integrat- 
ing out" the movable particles (up-spins), an effective interaction among the fixed particles 
(down-spins) is left. This, then, is a classical lattice gas with a complicated interaction. It 
resembles an Ising (not Heisenberg) model when U is large and it can be shown [KL] to 
have long range order at low temperatures (and no long range order at high temperatures) 
in two or more dimensions — as does the Ising model. 

C. The Surprising Hole. 

In the previous subsection we considered the half-filled band and showed a strong 
tendency to antiferromagnetism. The only important requirement on T was that it was 
real (i.e., no magnetic field acts on the orbital motion). Nagaoka [NY] made a surprising 
discovery about the case U = +oo, but with one hole (i.e., N — |A| — 1). Thouless [TD] had 
a similar result a bit earlier (cf. note 7 in [NY] and the discussion on p. 47 in [LE2]) but 
there is little doubt that Nagaoka's presentation of this particular result is clearer, more 
precise and applicable to certain non-bipartite lattices such as bee and fee. The Thouless 
approach uses a Perron-Frobenius argument that appears to be restricted to bipartite 
lattices; for such lattices it is true that fermions behave like bosons when U = -\-oo and 
there is one hole. Nagaoka, on the other hand seems to require a regular Bravais lattice, 
but this is not really necessary provided all t^y are nonpositive. The fully general result 
with a considerably simplified proof, was given by Tasaki [THl] . 

It is usually assumed in the condensed matter literature that the t^yS are nonnegative, 
but there seems to be neither a compelling reason for this assumption nor many examples 
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in which it can be verified (S. Kivelson, private communication). (For a bipartite lattice, 
however, one can have either sign with the help of the unitary operator (—1)^ applied 
to both spin | and |.) Notice that txy > puts the lowest kinetic energy at one point, 
namely /c = 0, in Fourier space. The condition txy < puts it at the edge of the Brillouin 
zone, and thus the lowest kinetic energy can be highly degenerate. The physical intuition 
is then quite diff'erent in the two cases — a fact that should not be lost sight of because of 
the hole-particle symmetry that holds for bipartite graphs. The proof in [THl] shows that 
negative, not positive txy is the natural mathematical assumption for this theorem. If this 
upsets anyone's physical proclivity, that is a pity. 

Theorem 2: If N = |A| — 1, ifU^ — +oo for all xeK, and ift^y < for all x, y, then 
the ground state has total spin S = N/2. This state is unique up to the trivial {N + l)-fold 
degeneracy if A satisfies a certain connectivity condition [THl]. 

The connectivity condition mentioned above is not stringent and it holds for all regular 
lattices in dimension greater than one (see also [AL]). Essentially it means that there are 
loops that permit nontrivial permutations of the particles. 

In the case that A is completely translation invariant, i.e., that all vertices of A are 
equivalent, as is the case on a hypercubic lattice on a torus, a different proof of the theorem, 
very similar to Nagaoka's, was given by Tian [TGI] and by Trugman [TS]. 

The obvious next question to ask is this: If there is more than one hole [N < |A| — 1) 
and U = +00, is the ground state totally ferromagnetic, i.e., is S = N/2? There can be 
no simple general theorem because numerical calculation on small systems show that the 
answer seems to be "no" and, at the same time, no simple pattern seems to emerge. Yet 
there are a few theoretical results, as follows. 

(i) The Nagaoka-Tian-Trugman method can be generalized for A = the 
D-dimensional hypercubic lattice with periodic boundary conditions and with t^y = t = 
nonpositive constant to show [TG2, TG3, TS, SQT] that the completely magnetized state 
energy, E{S = N/2; A), when compared to the actual ground state energy, -E(A), satisfies 



in the thermodynamic limit, A ^ oo, when the number of holes A^^ = |A| — A?^ is not too 
large. The best result is by Shen, Qiu and Tian [SQT], which gives (3.8) when Nh < |A|" 
with < a < 2/{D + 2). Note that we do not divide by |A| in (3.8), which thus truly 
represents the vanishing of a gap. The proof here is elegant and simple. However, one 
would expect (3.8) to hold as long as A^^/|A| — > as |A| — > 0. 

(n) If there are many holes, Nh/\A\ > 0.49 for the D = 2 or Nh/\A\ > 0.32 for the 
D = 3 hypercubic lattices, and A — > oo as in (i) then 



There really is an instability of the S = N/2 state with respect to one spin flip. This was 
proved by Shastry, Krishnamurthy and Anderson [SKA]; see also [SAl]. The estimate was 



lim E{S = N/2; A) - E{A) = 



(3.8) 



E{S = N/2; A) - E{S = N/2 - 1; A) 74 0. 



(3.9) 
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improved to Nh/\A\ > 0.29 by von der Linden and Edwards [LvE]; Hanisch and Miiller- 
Hartmann [HM] simplified the calculation (but not the estimate of 0.29). 

(iii) Several authors [DW, FRDS, SAl, TB] were able to prove, for a translation 
invariant hypercubic lattice model, that when there are two holes the ground state energies 
satisfy E{S = N/2 — 1; A) < E{S — N/2; A). However, assumptions have to be made about 
the relative lengths of the sides. Siito [SAl] extended this to 2, . . . , 6 holes for a bcc lattice. 

(iv) Siito [SA2] shows, as expected, that the energy splitting needed for demagneti- 
zation is, in any case, negligible. He proves that M(/5, /i, p), the magnetization per site in 
field h at density p = A^/|A| satisfies (in the thermodynamic limit) M(/?, h,p) tanh(/9/i) 
as p ^ 1. 

The results in (ii) and (iii) are achieved with a variational calculation. The value of 
E{S — N/2;A) is easy to find exactly because it equals the energy of spinless electrons, 
i.e., 

TV 

E{S = N/2;A) = J2^jiT) (3.10) 

i=i 

for any A and hopping matrix, T, and in which Ai(T) < A2(T) < ■ ■ ■ are the eigenvalues 
of T. Thus, the hard problem is to find a good variational function with S = N/2 — 1, and 
this appears to be extraordinarily difficult. Why? No one seems to know! And why is it 
so difficult to treat S = N/2 - 2? 

These results, (i)-(iii), show that one cannot expect S = A^/2 except when N = |A| — 1, 
but one can ask the following. 

Problem 2: With U = +oo, for which p := N/\A\ is it true that some (if there is 
more than one) ground state has S/\A\ > in the thermodynamic limit A — > oo? 

This brings us to two more open problems about the U = +oo case; the first is a 
corollary of the second. We take a large A and particles and suppose that the ther- 
modynamic limit A ^ oo with p = A^/|A| fixed is well defined. We set S = spin of the 
ground state (the maximum such spin if there is more than one ground state) and we set 
'S'max = A'"/2. We also assume that t^y > for all x and y. Then 

Problem 3: Prove or disprove that 

lim lim 5/5max = 0. (3.11) 

p— *-0 A^co 

Problem 4: Does there exist some number Pc> such that 

lim S/Sma^ = for all p < Pc? (3.12) 

A— i-oo 

The requirement that txy > is important. As we shall see in Sect. D, Miclke's work 
shows that there can be nice, periodic lattices in any dimension (such as the kagome lattice 
in two-dimensions) for which S = S'max for ah P < Pc with Pc > 0, thereby contradicting 
(3.11) and (3.12). To achieve this, however, one needs txy < 0. Perhaps (3.11) and (3.12) 
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hold in the case txy < if we replace "maximum such spin" by "average of such spins" , 
because Mielke's and Tasaki's examples have highly degenerate ground states with spins 
ranging from to N/2. 

Closely related in spirit to Theorem 2, but with an interesting, different proof, is the 
result of Chakravarty, Chayes and Kivelson [CCK]. They start with a half-filled band, 
N = |A|, and U large. Then they add or subtract one or two particles, and define En := 
E(N = |A| + n) for -2 < n < +2. They then define 

Ae = 2Ei — E2 — Eq 

Ah = 2E-i-E-2-Eo. (3.13) 

The interpretation of Ae is as a 2-particle binding energy, while that of A^ is as a 2-hole 
binding energy. The picture of Ag, for example, is that given two very large systems 
at half-filling, and given two extra electrons, is it energetically favorable to add the two 
electrons to one system (Ae > 0), or is it favorable to add one electron to each system 
(Ae < 0)? The former, Ae > 0, connotes pair binding. 

It is pointed out in [CCK] that Ae < and A/j < when U = 0, but they quote 
numerical studies showing that Ag > and A^ > for some U and some A. They prove, 
however, that in the limit U ^ 00, Ag < if all txy are nonnegative and A^ < if all 
txy are nonpositive. For a bipartite graph the sign does not matter (as long as all txy are 
positive or all are negative), and thus Ae < and Ah < in this case. They interpret 
this result to mean that the numerical positive binding results are only an intermediate U 
phenomenon but, strictly speaking, the (unlikely) alternatives Ag = and Ah = in their 
theorems would first have to be eliminated. 

D. Another Path to Ferromagnetism. 

The one-hole, U — +00 model is not the only one known to have saturated ferromag- 
netism (i.e., S = N/2). Mielke [Mil] and later Tasaki [TH2] and then both [MT] found 
interesting, but very special models with this property. 

The basic idea is to find a graph A and a hopping matrix T such that the lowest 
eigenvalue of T (call it Aq) is highly degenerate; in fact we want the degeneracy to be at 
least and, to be interesting, we want that to be of the order of |A|. Let A'^o denote the 
degeneracy of this lowest eigenvalue and denote the space of these eigenfunctions by Hq. 
Ux > for all a; G A, it is easy to see that if = A^o then 

(i) There is a ground state having S = No/2. 

(ii) The state is simply a determinant formed from the Nq vectors in Ho and its energy 
is A^oAo- 

(iii) The ground state is unique if certain additional conditions (known to be optimal) are 
met [Mil]. 

If A^ < A'o the ground state manifold will contain at least one state with S — N/2, 
but perhaps others as well (see [Mil]). 
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A comparison with Theorem 1 is interesting, but it is not clear whether or not it is 
misleading. Note that the S = \ \A\ — \B\\/2 result there for a bipartite graph and a half- 
filled band was somehow related to the | \ A\ — \B\ |-fold degeneracy of the zero-mode of T. 
Thus, a common feature is degeneracy, and it is often said that itinerant ferromagnetism 
is associated with atomic or kinetic energy degeneracy. But there are also important 
differences: 

(i) The spin in Theorem 1, while it might be proportional to |A|, is not N/2 = |A|/2. 

(ii) No fine tuning of T or of A was needed for Theorem 1. All that was needed was the 
bipartite structure, the positivity of the Ux^s and the reality of the t^^'s. 

Mielke's way of achieving the degeneracy A'^o is to start with some graph G and then 
to set A = L{G), the line graph of G, which is defined as follows. Make a mark in the 
center of every edge of G; those marks will be the sites of L{G). Two sites are connected by 
an edge in L{G) if the two edges of G on which they reside have a G-site in common. Note 
that L{G) is never bipartite, except for the trivial case of a ring. A well known example of 
a line graph is the kagome lattice. (Incidentally, kagome is not a person — it is a pattern 
of woven bamboo.) 

The hopping matrix T is defined to be — t < on every edge of A = L{G). Not only 
is there a restriction on the magnitude of t^y but we see, once again, that negative t^y is 
the natural sign — as it is for Theorem 2. 

Of course, as Mielke and Tasaki note, and as is also true in Theorem 2, one can convert 
positive txy into negative t^y by a hole-particle transformation on both spin | and spin J,. 
This does not alter the interaction or the total spin, S, but it changes N into N' = 2|A| —N. 
For positive, t^y, then, we can transform to the negative t^y situation and conclude that 
S = N'/2 when A^' < Nq. This translates into 5 = |A| - N/2 when N > 2|A| - Nq 
(and not <). This construction would yield saturated ferromagnetism, S — N/2, only at 
A = |A| (half-filled band), provided it were possible to achieve Nq > |A|; this is clearly 
impossible (unless T = 0), so the positive t^y choice does not yield the desired saturated 
ferromagnetism. However, it is still possible to have unsaturated ferromagnetism, i.e. 
1 > 2S/N 7^ in the thermodynamic limit. By taking N = 2|A| - Nq and 5' = |A| - N/2 
we have 2S/N = Ao/(2|A| —Nq). (One can reasonably argue, however, that the condition 
5 = |A| — A/2 is indeed saturated ferromagnetism because this value of the spin is the 
maximum possible one when N > |A|, given the constraints of the sytem and given the 
Pauli principle. The system is constrained by allowing only a limited number of states for 
the electrons, i.e., two per site. I leave the semantic resolution to the reader.) 

The homology of G determines Hq in a simple way. Pick any closed, self avoiding 
path in G of even length. This path corresponds, in an obvious way, to a closed path, P, 
of the same length in L{G). If, now, we take the vector 4'{x) = for x not a vertex of P 
and (p{x) = alternately -f-1 and —1 as we traverse the vertices of P, we see at once that 
(p is an eigenvector with eigenvalue —2t. It is also not hard to prove that —2t cannot be 
improved, i.e., Aq = —2t. Moreover, this construction yields all of Hq. 



14 



13/Nov/93 



In [TH2, MT] essentially the same result as Mielke's (with similar requirements on T) 
is achieved, but with certain decorated lattices with next nearest neighbor hoppings. The 
ferromagnetism of [TH2, MT] is shown to be stable under small change of the electron 
density, by using a "grand canonical ensemble" with a fixed electron density. It is also 
proved that there is a transition to a paramagnetic phase as one decreases the electron 
density. Thus (3.11) and (3.12) (with S replaced by its grand canonical average) are proved 
for these special graphs. 

In all cases, the eigenvectors in Ho can all be taken to have compact support, i.e., 
each (f){x) vanishes except on a finite set of sites of A, and each such set is independent of 
A once A is large enough to include the set. This property leads to the result [MT] that 
the effective Hamiltonian, in a suitable subspace of states that includes the ground states, 
can be written exactly as a Heisenberg Hamiltonian. 

4. The Flux Phase Problem 

Very little seems to be known rigorously about the effect of a magnetic field on the 
orbital motion, that is, if we set 



for some vector potential A. The spin of the ground state might well change. The known 
proofs of Theorems 1 and 2 fail when A ^ 0. If, indeed, the spin of the ground state 
changes then we have a new kind of "magnetic field — spin" interaction, brought about by 
the Pauli exclusion principle. Indeed, something similar is discussed in the one-dimensional 
context in [FKl]. 

The energy certainly does change, and it is by now well known that when U = the 
zero flux state does not give the lowest ground state energy. When A is a square lattice 
and when A?" = |A|, it is conjectured that the maximum (!) possible flux, namely tt in each 
square (^12^23^34^41 = ~1 around a square), is optimum. (For a survey of this question 
and for some rigorous results about it see [LL2].) The same conjecture has been made for 



Problem 5: Solve the flux phase problem for U ^ Q (or even for U — 0) for a 
half-filled band on a square lattice. 

The moral of this story is that the Pauli exclusion principle can really upset our ideas 
about diamagnetism. For one solitary electron, the imposition of a magnetic field raises 
the ground state energy. For many electrons it can and does lower the energy. It is trivially 
true, for example, that when U = +00 and N = |A|, or ?7 = anything and N = 2|A|, a 
magnetic field has absolutely no effect on the energy. 





(4.1) 
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5. Uniform Density Theorem 

The hole-particle transformation has remarkable consequences for the Hubbard model 
on a bipartite graph with a half-filled band, N = |A|. The ideas given below are well 
known to chemists — less so to physicists. They go back to Coulson and Rushbrooke [CR] 
for the U = case, and to MacLachlan [MA] who generalized them to many interacting 
models — including the Hubbard model as a special case. A simplified proof together with 
an extension to models involving explicity spin-spin interactions (such as the t — J model), 
and to the Falicov-Kimball model, is given in [LLM]. 

The results apply equally to three cases: (a) The canonical Gibbs state with N = |A|; 
(b) The grand-canonical Gibbs state with zero chemical potential; (c) The ground state 
with N = |A| and which is defined, in case of degeneracy, to be the /? — > oo limit of the 
canonical Gibbs state. The one-body density matrix p^{x,y) is the expectation value of 
^xa^ya State in question. 

Theorem 4: For a half-filled band on a bipartite lattice, the one-body density matrix 
for each of the above states satisfies (for each o" =t or J, ). 

Pa{x,y) = \5xy if x,yeA or x,yeB. (5.1) 

If xeA and yeB, nothing simple can be said. Note that the theorem does not require 
T to be real, and is thus one of the few theorems that applies to complex t^y. However, a 
true, physical magnetic field would also act on the electron spins and thereby vitiate the 
hole-particle symmetry needed for the proof. 

The complex case is a bit subtle, for it uses more than just a hole-particle transfor- 
mation (call it W). It also utilizes the nonlinear antiunitary map J that maps a vector 
tjj, considered as a polynomial in the c^^^s applied to the vacuum, into the vector tj;* cor- 
responding to the polynomial with complex conjugate coefficients. While J is nonlinear, 
JKJ is linear when K is any linear operator, and TrJKJ = (TrK)*. The antiunitary 
Y = JW satisfies Y = JW = WJ (and hence Y^ = 1) and Yc^^Y = Wc^^W. Most im- 
portant is the invariance of the Hamiltonian, YHY = H, which replaces the hole particle 
invar iance, WHW — H, which fails for complex T. 

What is the significance of this result? It seems to contradict the conjecture in Problem 
1 that there can be antiferromagnetic, i.e., staggered, LRO. Such an ordering can occur 
only in the thermodynamic limit and it has the property that for every state with ordering 
titi ■ ■ ■ there is a state with ordering jfj,! • • •. The point about Theorem 4, applied to 
the thermodynamic limit, is that for every state with ordering there is an equally good 
state (obtained by changing boundary conditions) with the opposite ordering, and on the 
average each site will have the same density for each spin value. In other words, there is no 
way to adjust the potentials Ux or the hopping matrix txy in a clever way so as to enhance 
the occupation of certain sites — in a manner independent of boundary conditions. This 
stability is remarkable and, although it is not true for real materials, the theorem hints 
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at some kind of remnant stability that might transcend the overly idealized assumptions 
needed for its proof. 

Problem 6: Is there any residue of (5.1) when N ^ |A|? 
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